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We consider a model fluid with long-ranged, (dispersion) interparticle potentials confined 
between competing parallel walls. One wall is solvophilic and would be completely wet at bulk 
liquid-gas coexistence fi^o while the other is solvophobic and would be completely dry at fi = 
When the wall separation L is large and the system is below the bulk critical temperature Tc and 
close to bulk liquid-gas coexistence, a 'delocalized interface' or 'soft mode' phase forms with a liquid- 
gas interface near to the centre of the slit; this interacts with the walls via the power-law tails of the 
interparticle potentials. We use a coarse-grained effective Hamiltonian approach to derive explicit 
scaling expressions for the Gibbs adsorption, F, the surface tension, 7, the solvation force, fs and the 
total susceptibility x- These quantities depend on the dimensionless scaling variable (L/a)^ /SS/i, 
where j3 = (ksT)'^ , a is the diameter of the fluid particles and Sfj, — fi — fico is the chemical 
potential deviation from bulk coexistence. Using a non-local density functional theory (DFT) we 
calculate density profiles for the asymmetrically confined fluid at different chemical potentials and, 
for sufficiently large L, confirm the scaling predictions for the four thermodynamic quantities. Since 
the upper critical dimension for complete wetting with power-law potentials is < 3 we argue that 
our (mean-field) scaling predictions should remain valid in treatments that incorporate the effects 
of interfacial fluctuations. As the wall separation L is decreased at fico we predict a capillary 
evaporation transition from the 'delocalized interface' phase to a dilute gas state with just a thin 
adsorbed film of liquid-like density next to the solvophilic wall. This transition is connected closely 
to the first order pre-wetting transition which occurs at the solvophilic wall in the semi-inflnite 
system. We compare the phase diagram for the competing walls system with the phase diagrams 
for the fluid confined between identical solvophilic and identical solvophobic walls. Comparisons are 
also made with earlier studies of asymmetric confinement for systems with short-ranged forces. 



I. INTRODUCTION 

Understanding the phase behaviour and equilibrium 
structure of confined fluids is a subject of intrinsic inter- 
est to the statistical mechanics communities and is also 
a key ingredient in elucidating the microscopic dynamics 
and flow properties of fluids in microfluidic or nanoflu- 
idic devices. When a simple (atomic or molecular) fluid 
is confined in narrow pores, or channels or capillaries, 
where the size of the confining dimension is of the order 
of nanometers, the effects of both finite pore size and 
substrate fluid interactions (adsorption) can have a pro- 
found influence on phase behaviour giving rise to phe- 
nomena that have no direct counterpart in bulk [l|, 
For a real fluid conflned by a real structured substrate or 
in a real porous solid details of the substrate (wall) po- 
tentials matter. The nature of the wall-fluid potentials, 
the roughness of the walls, the atomic corrugations and 
the geometry of pores or grooves all play a role in deter- 
mining the microscopic structure of the confined fiuid. 
However, when it comes to phase behaviour one might 
hope that certain generic features would be captured by 
simple models of confinement. This is the motivation be- 
hind the majority of theoretical and computer simulation 
efforts to ascertain the nature of phase transitions of fiu- 
ids confined in idealized geometries; these are often pla- 
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nar (slit-like) or cylindrical model pores with structure- 
less walls representing the substrate [l|, . Indeed much 
of the progress in the field has resulted from detailed 
studies of even simpler models namely a lattice gas or a 
nearest-neighbour Ising film of flnite thickness D where 
the spins at opposite walls are subject to (local) surface 
flelds Hi and Hu in addition to the external, bulk mag- 
netic fleld H. Phenomena such as capillary condensation 
and evaporation, prewetting and the interface delocal- 
ization transition have been investigated and elucidated 
in Ising model studies Q following on from the work of 
Fisher and Nakanishi in the early eighties [1, HI . In re- 
cent years some focus has shifted to soft matter systems, 
e.g., colloids and polymers, where the mesoscopic length 
scales of the 'particles' are such that the details of the 
atomistic structure and roughness of the conflning walls 
should not be relevant. Moreover the large size of the 
colloidal particles means that these can often be tracked 
in real space and in time using confocal and video mi- 
croscopy thereby revealing information about structure. 
Of course, theoretical techniques developed for atomic 
fluids are easily adaptable to colloidal systems — the latter 
having the advantage that by applying suitable chemistry 
the effective fluid-fluid and wall-fluid interactions can be 
tailored. For an informative overview of models of soft 
matter, in particular colloid-polymer mixtures, conflned 
in thin film geometry, i.e., between two parallel walls, 
see Ref. Q. The phase behaviour of symmetrical poly- 
mer blends in the same confining geometry is reviewed 
in Ref. [3]. 
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In the present paper we consider confinement of a sim- 
ple fluid between 'competing' planar walls one of which 
is solvophobic — it is completely wet by the gas phase — 
whereas the other wall is solvophilic and is completely 
wet by the liquid phase. The phase behaviour in such a 
system is very different from that when the two walls are 
identical 

Parry and Evans [1, 0] drew attention to the novel 
phase equilibria that may arise when a fluid (or Ising 
magnet) is confined between two perfectly antisymmet- 
ric walls. For temperatures below the bulk critical tem- 
perature (T < Tc) they identified a single 'soft mode' 
phase characterised by a liquid-gas (-1 — spins in the 
Ising case) interface which lies parallel to and at the mid- 
point between the walls when the fluid in the reservoir 
is at bulk liquid-gas coexistence. In the Landau the- 
ory studied by Parry and Evans [1, the forces are 
short-ranged and the liquid-gas interface fluctuates freely 
about the centre of the slit with an exponentially diverg- 
ing transverse correlation length ^|| ~ e'*^/'* and the to- 
tal susceptibility diverges as e"^^/^, where L is the wall 
separation. In the original, mean-field Landau analysis 
the length scale is the bulk correlation length: 
= ^b- Subsequently Parry and co-workers |10l - [l3 | 
used a generalized effective Hamiltonian approach to ar- 
gue that ~ ^b(l+ti;/2) where the wetting parameter 
Lo is defined as = ksT j (Air^'^^ig), ^ig being the in- 
terfacial tension for the free liquid-gas interface or the 
interfacial stiffness for the Ising case and ^f, is the true 
bulk correlation length. This 'soft mode' phase persists 
until a critical temperature Tcl where Tcl < Tw, the 
critical wetting transition temperature at a single wall 
(for perfectly antisymmetric walls the single wall drying 
transition temperature To is equal to Tw)- Below T^l 
symmetry breaking occurs and there is two phase coex- 
istence between states where the liquid-gas interface is 
bound to either one of the walls. The transition at Tcl is 
sometimes termed an interface localisation-delocalisation 
transition and can be first or second-order depending on 
the single wall wetting and drying transitions. 

This unusual behaviour for a system subjected to an- 
tisymmetric surface fields has been studied in several pa- 
pers by Binder et. al. [l3l - [l5| using Monte Carlo simu- 
lations of a nearest-neighbour lattice gas model (equiv- 
alent to an Ising model in magnetic terminology). On 
decreasing the temperature T below the bulk critical 
temperature Tc the layer magnetisation m„ (or density 
Pn = [1 — Tnn)/2]) profiles (see, e.g.. Fig. 1 in Ref. [T^ ) 
show the gradual formation of a -I — interface between 
the two bulk phases. The presence of this interface is also 
evident in the layer susceptibility profile Xn — dnin/dH, 
where H is the external magnetic field equivalent to the 
reservoir chemical potential deviation dp = fi — pco in 
the lattice gas model. Xn develops a pronounced peak 
in the centre of the film. Note that in the Ising simula- 
tions the wall separation is denoted D rather than L, and 
the surface fields, acting only on layers 1 and D, satisfy 
Hn = —Hi. For large L the susceptibility is large at 



H = because the interface is almost freely floating — 
the free energy cost for shifting the whole interface is 
exponentially small. The simulations of Binder et. al. 
confirmed the Parry and Evans prediction [9] of a suscep- 
tibility which diverges exponentially with L. Binder et. 
al. |13l - [l5| also studied the phase transition at Tcl , where 
the symmetry is broken and the interface becomes bound 
to one of the walls rather than delocalised in the centre of 
the slit. For critical wetting, in agreement with Parry and 
Evans [1, Q , they found that Tcl lies close to but below 
the single wall wetting transition temperature T\y. As 
expected this interface localization-delocalization transi- 
tion has two-dimensional Ising character. Monte Carlo 
simulation results for the Ising model are reviewed in Q 
where they are compared to mean field predictions. 

Following these developments in the early and mid 
nineties the phase behaviour of other types of fluid 
confined between competitive walls was investigated. 
Miiller and co-workers studied binary (symmetric) poly- 
mer blends subject to asymmetric surface fields using 
Monte Carlo and self-consistent field approaches — see 
e^., Ref. [l6l |. A summary of such studies is given in Ref. 
[31 ■ It is well-known that sterically stabilized colloidal 
suspensions with added non-adsorbing polymer exhibit 
fluid-fluid phase separation into phases that are rich or 
poor in colloid. The bulk properties of such mixtures 
are often described by the Asakura-Oosawa-Vrij (AO) 
model. De Virgiliis and co-workers have investigated the 
phase behaviour of the AO model fluid between asym- 
metric conflning walls where one wall is completely wet 
by the colloid rich phase and the other by the polymer 
rich phase 

For completeness we also mention recent studies of 
the (confined) Ising strip in two dimensions subject to 
various surface fields [19-23'], including the antisymmet- 
ric Hd — —Hi case. These reveal 'pseudo' transitions 
in the reduced dimensionality that are analogous to the 
real transitions occurring in three dimensions. There are 
also recent Monte Carlo studies of Ising films subject to 
general asymmetric (local) surface fields [13] • Note that 
in Ref. [22\ simulations were performed for the nearest- 
neighbour Ising strip subject to l ong -ranged anti- 
symmetric surface fields while Ref. 19] performed density 
matrix renormalization group calculations for the same 
strip subject to various power-law antisymmetric surface 
fields. Ref. [2^ presents results for the 'soft mode' phase, 
discussing in detail the effects of capillary wave fluctua- 
tions. These effects are much more pronounced in spatial 
dimension d = 2 than in d = 3. 

Here we study the behaviour of a simple fluid subject 
to competing walls at a single temperature T in the 'de- 
localised interface' or 'soft mode' phase, i.e., Tcl < T < 
Tc- In contrast to previous studies of fluids [1, [13, [11], 
where the interaction potentials were short-ranged, we in- 
clude the full tail, —r~^, of the Lennard- Jones 12-6 pair 
potential in the fiuid-fluid and wall-fiuid interactions. As 
a result the centrally located liquid-gas interface inter- 
acts with the two walls via the power-law tails of the 
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interparticle potentials. The overall phase behaviour is 
expected to be similar to that described above for short- 
ranged forces but with quantities such as the suscepti- 
bility X ^"^d transverse correlation length ^|| diverging 
algebraically rather than exponentially 

In Section|ll]we employ an effective potential approach 
in which the excess grand potential is written as a func- 
tion of the thickness of the drying film I (this is equiva- 
lent to the 'slab approximation' of Parry and Evans, see 
Section 2.3 in Ref. [ol]). Our findings apply to the gen- 
eral case of competing walls, where one is drying and the 
other wetting but the walls are not necessarily perfectly 
antisymmetric. We consider state points both on and 
away from bulk coexistence ^^o- In Section III Bl expres- 
sions are derived for various thermodynamic quantities 
including the Gibbs adsorption, F, surface tension, 7, sol- 
vation force, fs and the total susceptibility x- For large L 
these are given by scaling functions of the dimensionless 
product {L / a)^ /36 ^ where a is the hard sphere diameter 
of the fiuid particles, 6 fi is the chemical potential devia- 
tion from bulk liquid-gas coexistence and f3 — [kBT)~'^. 
At bulk coexistence, 5^ — 0, the total susceptibility is 
predicted to diverge as x i"*- The solvation force, i.e., 
the excess pressure due to confining the fluid, is predicted 
to be repulsive: fs decreases with wall separation as , 
for large L. 

Capillary condensation (or evaporation) is a well 
known phenomenon for a fluid confined between two 
identical parallel walls. In Section fll CI we display phase 
diagrams (at constant temperature T = O.STf) for dif- 
ferent systems obtained using the effective potential ap- 
proach. First we consider two identical parallel solvo- 
phobic walls and determine the line of evaporation tran- 
sitions as a function of (5/i and inverse wall separation 
. In the second system the walls are both solvophilic. 
A single wall {L — 00) undergoes a first order wetting 
transition at Tw < 0.8Tc and the accompanying pre wet- 
ting transition can also occur in the confined fiuid where 
it competes with capillary condensation. We show that 
both capillary condensation and evaporation are possible 
in a competing walls system. For the particular system 
which we study in Sec. lIVI using density functional theory 
(DFT) the effective potential approach predicts a capil- 
lary evaporation transition which occurs for increasing 
6fi as the wall separation L is decreased. The line of 
evaporation transitions crosses bulk coexistence /i = fico 
at L = L'^yap- 

In Section Hill we describe sum rules that relate the sol- 
vation force to the contact densities at the walls. These 
are used in checking the accuracy of the DFT calcula- 
tions. 

Section IIVI describes numerical results obtained using 
a non-local DFT. The wall-fiuid potentials have been 
chosen so that at T = O.STq the two walls are 'anti- 
symmetric', i.e., the Hamaker constant [defined in Eq. 
(|B6|) ]. for wetting at the solvophilic wall [62 in Eq. ^] is 
equal to that for drying at the solvophobic wall [61 in Eq. 
(O]. Results are given for four different wall separations 



L/cr = 50, 100, 250 and 500 at temperature T = O.STc- 
We find good agreement with the effective potential pre- 
dictions of Section ini i.e., scaling is obeyed for all the 
thermodynamic quantities we consider, provided L is suf- 
ficiently large. Density profiles of the two coexisting 
states at L'^°ap — 11.4cr, obtained from DFT, are also 
presented. At L = lOOcr capillary evaporation occurs on 
the gas side of bulk coexistence at (3Sfj, = —4.13 x 10"'^ 
and we display the coexisting density profiles at this tran- 
sition. These are almost identical to the profiles at the 
prewetting transition. 

Our results are discussed in Section |Vl where we also 
speculate on the phase diagram as temperature is varied. 
We relate our findings to simulation and theoretical stud- 
ies of other model fluids and comment on their relevance 
for experimental results for confined water. 



II. EFFECTIVE INTERFACIAL POTENTIAL 
DESCRIPTION AND ITS PREDICTIONS 

Following from the pioneering work of Frumkin and 
Derjaguin in the 1930's (ii isig htfuUy reviewed in recent 
articles by Henderson [2^, [2^]), theories of wetting and 
condensation phenomena often express the excess grand 
potential of a confined fiuid as a function of a single varia- 
tional parameter, I, the thickness of a wetting film. Here 
we extend and quantify this effective interfacial poten- 
tial approach, introduced in Ref. [9| for asymmetric con- 
finement, specializing to the case of power- law attractive 
potentials. 



A. The model 

The specific model fluid we consider here has a pair 
potential consisting of a hard core of diameter a plus an 
attractive tail, taken to be the attractive part of the full 
(12,6) Lennard- Jones (LJ) potential: 



(t>att{r) 



4e 



r , 



(1) 



where r^m — 26cr and e is the LJ well-depth. We 
shall describe this fluid in the context of a density func- 
tional theory (DFT) — see Section IIVI — that treats the 
hardcore repulsion within fundamental measures theory 
and the attractive interactions ([T]) in mean-fleld fashion. 
In this DFT approach the bulk critical temperature is 
ksTc/e = 1.415. The fluid is conflned between parallel 
walls separated by a distance L, see Fig. [TJ The wall- 
fluid potentials are full Lennard- Jones (9,3), obtained 
by integrating the individual Lennard- Jones (12, 6) wall 
particle-fluid particle interactions with the same form 
as (H)) over the volume of the wall, see Eqs. (|A1I) - (|A2I) . 
There is hard repulsion at the contact surface of the wall 
with the fluid. The strength of the potential exerted by a 
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FIG. 1. Diagram to show the sharp-kink approximation for 
the density profile p{z) of the fluid in an asymmetric slit. 
The walls are comprised of uniform-density blocks of particles, 
separated by a distance L. The wall-particles interact with 
the fluid particles with the same potential as the fluid-fluid 
interactions. A fluid particle is shown at the plane of contact 
of the fluid with walli, i.e., at 2: = 0. A wall particle (shaded) 
is shown at the furthest extent of walli towards the fluid. 
There are excluded volumes, of width dw = a, in which the 
density is zero between the fluid and the regions of constant 
wall density p^^ and pw2- In the sharp-kink approximation 
there is a film of gas, width I, with the coexisting bulk gas 
density pg next to (drying) walli; pmi < Pg- The remainder 
of the slit is filled with a film of liquid, width L ~ I, with the 
coexisting liquid density pi. Wall2 is wetting; p^^ > Pi- The 
density profile is discontinuous at the gas-liquid interface at 
z = l. 



wall is varied by changing the density of wall particles p^, . 
The wall potential is made repulsive by choosing a nega- 
tive value for p^, . For simplicity the range of the wall- fluid 
interparticle potential has been chosen to be the same as 
that of the fluid-fluid potential, i.e., a^f — cr. There is 
an excluded volume of width dw (in this case dw — a), 
measured between the centres of the particles in the sur- 
face of the wall and the centres of the fluid particles at 
the contact surface of the fluid density profile [131, as il- 
lustrated in Fig. [TJ The total planar external potential 
on a fluid particle confined between the two parallel walls 
is V{z) = Vwi{z) + Vw2{L — z). Vw^{z) is the potential 
exerted on a fiuid particle at z by a single walli, with 
the walli-fluid contact surface at z = 0. V^^iL — z) is 
the potential due to the second wall (wall2) which has its 
contact surface with the fluid aX z = L. The wall den- 
sity piuj can be chosen so that in the semi-inflnite system 
consisting of walli in contact with the liquid the wall is 
completely dry (wet by the gas) at liquid-gas coexistence 
and conversely the density of wall2, Pw2^ can be chosen to 
ensure complete wetting by liquid. The fluid between the 
two walls is in contact with a reservoir with a chemical 
potential p. and temperature T, i.e., the system is in the 
grand canonical ensemble. 



B. One wall wet and the other dry 

We focus first on the asymmetric confinement where 
walli is completely dry and wall2 is completely wet in 
the limit L — > 00. In the sharp-kink approximation the 
fluid density is constant on both sides of a step function 
situated at a distance I from walli, as shown in Fig. [T] 
Using this approximation we can write down the excess 
grand potential per unit area A for the system at bulk 
liquid-gas coexistence p, — p-co- 

A " A 

+ Wyj^g/gl{l) + W^^l/lg{L ~ I) 

+Wwig/w2l{L), (2) 

where p = p{fJ-, T) is the pressure in the reservoir, V is 
the accessible volume and I is the thickness of the gas 
film. The first three terms correspond to the individual 
surface tensions of the three interfaces: walli-gas, wall2- 
liquid and gas-liquid. The next terms give the interaction 
potentials between pairs of interfaces. Using the sharp- 
kink approximation (see Appendix [B]) we find: 

tW^^ig/gKO = + (3) 

"^-'/'«(^-^)= (L^+(L^' 

Wwig/w2l{L) = ^ + ^, (5) 

where the temperature dependent (Hamaker) coefficients 
are 



hi 


= {pg - Pw,){pi - Pg)nea^/3, 


(6) 


b2 


= {Pi - Pw2){Pg - p^T^ecrVS, 


(7) 


bs 


= {PwiPi + Pu)2Ps)7i'eo-*^/3, 


(8) 


Cl 


= 2apwi{pi - Pg)'Kea'^ /i, 


(9) 


C2 


= 2apnj2{Pg - Pi)T^<^<y^ /i, 


(10) 


C3 


^ ~2a{p.u,iPi + p.u,2Pg)TTe(j^ /i. 


(11) 



Pg and pi are the densities of the coexisting gas and liquid 
at temperature T. For simplicity we have set the wall- 
fiuid parameters twif = ^W2f = e in the above equations 
but since it is the products Pwi^wif and Pw2^w2f which 
determine the strength of the wall potentials the expres- 
sions are still completely general. The lowest order coef- 
ficients, &i, 62 and 63 are expected to be correct beyond 
the sharp-kink approximation but we anticipate that ci , 
C2 and C3 will be affected by the microscopic details of the 
density profile. The direct interaction between the two 
walls [ not due to the presence of the fiuid) is not included 
in the excess grand potential As already implied, in 
the analysis which follows the temperature and wall den- 
sities are chosen such that bi and 62 are both positive, 
i.e., in the semi- infinite systems consisting of a wall in 
contact with the liquid at bulk coexistence, walli would 
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be completely dry and conversely wall2 would be com- 
pletely wet. 

If the system is not precisely at bulk liquid-gas co- 
existence then there is an additional term in the excess 
grand potential because one of the (bulk) phases present 
is metastable: 



Vtex{l\L,n) = Vtex{l\L,llco) + {p - Pm)Vm 



(12) 



where Pm and Vm are the pressure and volume of the 
metastable phase at the same chemical potential. On the 
liquid side of coexistence, i.e, ^ — ^co = > 0, and for 
small 5^ the pressure difference is (p — pm) = <5/z(p; — Pg) 
and the volume of the metastable gas is Vm = I A. li dfi 
is negative then it is the liquid phase which is metastable 
and (p - Pm)Vm = \Sn\ipi - Pg){L - l)A. 



1. Film thickness and adsorption 



that in the perfectly antisymmetric case, where leq = L/2 
for (5/i — 0, A[0] = 1/2 which requires bi = 62 at the tem- 
perature of interest. 

We define the (dimensionless) Gibbs excess adsorption, 
r, in terms of the equilibrium density profile p{z): 



r(L,/.) 

A 



0+ 



dz[p{z) - p], 



(17) 



where p = p{ii/T) is the reservoir density. On the liquid 
side of bulk coexistence T is negative and the gas film 
thickness in the sharp-kink approximation is given by 



A{PI - Pg) 



Sfi > 0. 



(18) 



On the gas side of bulk coexistence F is positive and the 
excess adsorption is proportional to the thickness of the 
liquid film, L — l^q, 



Substituting Eqs. ([2])-® into Eq. (HU we obtain the 
excess grand potential for the system on the liquid side 
of coexistence. 



lwig{P-co) + 7u.2i(Mco) + llg 



Cl 



{L-iy L2 

C2 , C3 



+ 7^ +<5M(Pi - Pg)l 



P {L~lf L3 
+0 \l-\{L-l)-\L-~^] 



(13) 



The smoothness of the actual profile is responsible for 
the leading order corrections, etc., to this effective 
potential; recall ci, C2 and C3 are the sharp-kink values 
for the coefficients. The corresponding equation for the 
gas side of coexistence is obtained by replacing the term 
(5/i(/9; — pg)l with 5^{pi — Pg){l — L). The equilibrium gas 
film thickness leg, is found by minimising Eq. (jl3|) w.r.t. 
I, 



2bi 



2&2 



3ci 



3C2 



4 +Sfl{pi-Pg) = 0. 

(14) 

Higher order terms are not displayed. 

In the limit L — > cx), ?eg — >■ 00 we can neglect terms 
of order and higher. Then the equilibrium gas film 
thickness can be expressed in terms of a dimensionless 
scaling function A: 



leq{L,p) 



= A [|3S^,{L/af] , 



(15) 



where A satisfies 
26i 2&2 



A3 (1-A)3 



6pL^ipi-Pg), 6p>0. (16) 



With 6fi replaced by \6p\ the same scaling function ap- 
plies on the gas side of coexistence, i.e., for 6p < 0. Note 



A[Pl - Pg) 

Although the excess adsorption as defined in Eq. (|T7)) is 
discontinuous at Sp = this is because of the jump in the 
reservoir density. In the sharp-kink approximation we set 
p = Pg for dp < and p = pi for 5p > 0. The profile 
and the thickness of the drying film l^q are continuous at 
5p = 0. 

That the adsorption, and other thermodynamic quan- 
tities to be described below, should be functions of the 
scaling variable f3Sp{L/(T)^ is consistent with the heuris- 
tic scaling ansatz introduced by Parry and Evans, see 
Sec. 4 in Q, in the particular case of dispersion forces. 



2. Surface Tension 

Substituting the equilibrium gas film thickness into the 
excess grand potential (|13p . we find that the equilibrium 
surface excess free energy per unit area, i.e., the surface 
tension, depends on a scaling function, S: 



l{L,p) 



where 



^ex ij-eq] L, p) 

^ A 

= IwigiPco) + lw2l{Pco) + llg 

+ -^S [p5p{L/cjf] 



(20) 



(21) 



(1-A)2 
+6pipi - Pg)AL\ (22) 



when 5p > Q and 



E [p5p{L/af] = + 7T-^ + &3 



A2 (1-A)2 
+ \Sp\{pi-pg){l~K)L\ (23) 
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when Sii < 0. The scaling function E has dimensions of 
energy. It is straight forward to show that (|2ip and ([2 
satisfy the Gibbs adsorption sum rule 



r_ f dj{L,^i) 



(24) 



which also diverges as L^. Away from the region near 
to the liquid-gas interface the density profile is expected 
to be largely unaffected by small changes in the chemi- 
cal potential. Correspondingly the local susceptibility is 
expected to be smaller in these regions than for z ~ leq- 



3. Susceptibility 

The total susceptibility x{L,fJ-) measures the response 
of the confined fluid to changes in the chemical potential 
fi. It is defined as 



Within the sharp-kink approximation 



XiL,^i) = (pi - Pg) { 



(25) 



(26) 



which is continuous at 6p ~ 0. From Eqs. P31 [T51 [TO)) 
we find to leading order 



xiL,p) = -ip,-p,rLU-± 



2 .4 / 65i , 662 



A4 (1-A)'= 



(27) 



= C [/36p{L/af] 



where C \^/3Sp,{L/a)^'\ is a scaling function with dimen- 
sion length"^. energy" The susceptibility at bulk coex- 
istence (5/i = is 



x{L,Pco) 



-L\PI-P,? 



In the perfectly antisymmetric situation, bi = 62 

-L\pi-Pg? 



(28) 



xiL,p,co) = 



1926i 



(29) 



The local susceptibility measures the change in the 
fluid density at each point in the density profile p(z) as 
the chemical potential p is varied: 



x{z;L,p) 



f dpjz) 
\ dp 



(30) 



The local susceptibility in the region of the liquid-gas in- 
terface is expected to be large in magnitude as the liquid- 
gas interface shifts with Sp whilst maintaining its shape. 
An increase in the chemical potential p causes a decrease 
in the gas film thickness, l^q — > l^q + 6leq, and the den- 
sity profile near the liquid-gas interface is translated by 
a (negative) amount Sleq in the z-direction. The local 
susceptibility near to the liquid-gas interface is therefore 



xiz;L,p)^[^ 



dpjz) 
dz 



(31) 



^. Solvation force 

An important quantity in the theory of confined fiuids 
is the solvation force, fs{L,p), defined as the force per 
unit area, arising from the presence of the fiuid, that must 
be exerted on the walls to maintain them at separation 
L. fs does not include any direct interaction between the 
walls. Although fs is referred to as a "force" it has di- 
mensions of pressure and thermodynamically should be 
thought of as the excess pressure in the fluid due to con- 
finement, i.e.. 



fs{L,p)^-j 



d^{leq]L,p) 

dL 



p- 



(32) 



It can be calculated from the change in surface tension 
j{L,p) with wall separation L, 



fs{L,p) 



dl{L,p) 



dL 



(33) 



ti,T,A 



The solvation force between our asymmetric walls is 
found from (fT3l)and (HSJ; to leading order 



fs{L,p)^j^Fs[l3Sp{L/af] 

where Fg is another scaling function given by: 
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Fs [/36p{L/ar] = - 



2b. 



(34) 



(35) 



Fg has the dimensions of energy 



C. Capillary condensation, evaporation and 
prewetting 

Before we examine possible phase transitions in the 
fluid confined between asymmetric walls it is helpful to 
recall the transitions that can occur when the two walls 
are identical. Firstly we consider the capillary evapora- 
tion transition that may take place in a fiuid confined 
between two solvophobic walls. For a fluid with 6p > 
confined between two planar drying walls we equate the 
excess grand potential of the evaporated state filled by 
gas, rif^(L, p), with that of the state in which the slit is 
filled mainly with liquid apart from a layer of gas, thick- 
ness leq at each wall, nl,^{leq; L, p). Within the sharp- 
kink approximation the excess grand potentials of these 
two states are given by: 



^UL,P) 
A 



2jwtg{p) + 5p[pi - pg)L + 0{L' 



(36) 
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FIG. 2. Phase diagram for a fluid confined between two iden- 
tical solvophobic walls at T = O.STc- The wall-fluid inter- 
particle potentials are given by Eq. (|A2[) with pwHwif = 
— l.lOSeo-""^. This choice of parameters ensures that 6i > 
and therefore complete drying would occur at an isolated 
walli. Eq. (|38|l is used to calculate the coexistence line be- 
tween the two phases. The sharp-kink value was assumed for 
the coefficient bi and the liquid-gas surface tension yig was 
calculated using DPT. 



and 



"^IwigifJ-) + "^llg +'^+ "^Sflipi - pg)l 



0[{L~2iy 



0{L- 



(37) 



Again the equilibrium thickness leq is obtained by min- 
imising w.r.t. I. jwigifJ-) is the surface tension of the 
wall-gas interface and bi is the Hamaker constant for dry- 
ing at a single wall, given by Eq. (jH]). The corrections of 
order in Eqs. (p6|) and ([37| refer to the interaction 
potential between the two wall-gas interfaces; those of 
order (L — 2/)^^ in (|37l) arise from the interaction po- 
tential between the two gas-liquid interfaces. Capillary 
evaporation occurs when the two excess grand potentials 
are equal: 



evap 



3leq (Sp-)] 



(38) 



where L^yap is the wall separation at evaporation. The 
factor of 3 in the denominator reflects directly the effects 
of power-law (dispersion) forces — see [1^. For short- 
range fluid- fluid and wall- fluid potentials {pi — Pg)S^ = 
2jig/ [Levap — 2Zeq((5/i)]. Figure [2] displays the phase di- 
agram, calculated using Eq. ([55]) . for the fluid between 
identical, parallel solvophobic walls as the wall separation 
L and chemical potential deviation from coexistence dp 
are varied. As dp increases the wall separation at evapo- 
ration decreases reflecting the increased free energy cost 
(oc dp) of a volume of fluid at the gas density. 

More complex phase behaviour is predicted for our 
fluid conflned between two identical solvophilic walls with 
the same wall-fluid potential as wall2 in the asymmetric 
system. In order to understand the transitions that oc- 
cur in the confined fluid it is instructive to first consider 



the wetting transition that would take place in the semi- 
infinite system, L = oo, at wall2 on increasing the tem- 
perature along the bulk coexistence curve on the gas side 
p — p~g(T). In contrast to the solvophobic walli which 
is dry at all temperatures, at wall2 there is a first-order 
wetting transition at a temperature Tiy, i.e., the thick- 
ness of the adsorbed liquid film jumps from being flnite 
below T\Y to infinite (macroscopic) above Tyy- Above 
Tw the first-order thin-thick transition occurs off coex- 
istence at a chemical potential deviation from coexistence 
Sppw{T) = Ppw(T) — Pco < 0. In this 'prewetting' transi- 
tion the film thickness jumps from a thin to a thick but 
finite value. Prewetting persists up to the prewetting 
critical temperature T^^ . The temperature at which we 
perform our DFT investigations, T = O.STc, is above 
the wetting temperature but below the prewetting criti- 
cal temperature, i.e., Tw < O.STc < . Consequently 
we observe a thin-thick film transition in the semi-infinite 
system. This prewetting transition may also take place in 
the symmetrically confined fluid, at a very similar value 
of chemical potential to that of the transition at a single 
wall when L is large. The excess grand potential of the 
thin film state is: 



^lr{L,p) 

A 



2j,y,g{p) + 0{L-') 



(39) 



where ^w^gip) is the surface tension of the wall2-gas in- 
terface (the thin liquid film state) for the semi-inflnite 
fluid at a single wall. The excess grand potential of the 
thick liquid fllm state is: 
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thick 

v-^'-'-^-27^,,(m) + 27,, + 72 



{luq;L,p) 



A 



2&2 



(40) 



+2\5p\{pi - pg)luq + 0[{L - 2luqy 



where 7^)2/ (/i) is the surface tension of the wall-liquid in- 
terface, liiq is the thickness of the wetting (liquid) fllm 
and &2 is the Hamaker constant for wetting at a single 
wall, given by Eq. The thin-thick film transition oc- 
curs when rtl'l™{L,p) = Q.^^J-'^^{liiq;L,p). However the 
thin-thick fllm transition is now in competition with cap- 
illary condensation. The excess grand potential of the 
condensed liquid state is: 



ntconc 
'''•ex 



{L,p) 



A 



2^^,i{p) + \Sp\{p,- pg)L + 0{L-^). (41) 



The equilibrium state at any given wall separation L and 
chemical potential p is the state with the lowest excess 
grand potential. The phase diagram, calculated using 
Eqs. dSSl), (HD) and (gT]), is shown in Fig. H At con- 
stant (large) L, increasing the chemical potential results 
in two first-order transitions on the approach to bulk 
coexistence — first a thin-thick fllm prewetting transition 
at Sp Sppw, independent of L in the present approx- 
imation, and then a transition from a state with thick 
wetting fllms at both walls and gas in the centre to a 
condensed state in which the slit is completely fllled with 
liquid. At small L, i.e., a/L > 0.00465, there is a single 
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FIG. 3. Phase diagram for a fluid confined between two identi- 
cal solvophilic walls at T = O.STc- The wall-fluid interparticle 
potentials are given by Eq. HA2|) with pn>2^w2f ~ f.723eCT~''. 
This choice of parameters ensures that 62 > and therefore 
complete wetting would occur at an isolated wall2. The co- 
existence lines between the difi^erent states were calculated 
using Eqs. (|39p . (|4ip and pif) . assuming the sharp-kink value 
for the coefficient 62- DFT was used to calculate the sur- 
face tensions for the liquid-gas 7;g, wall2-liquid 71031(^00) and 
wall2-gas "iw2g{l^(^o) interfaces. Note that the wall2-gas surface 
tension 7^39(^00) is the excess grand potential per unit area 
of the thin film state and refers to a metastable state at bulk 
coexistence fj,co- The coexistence values for the wall2 -fluid 
surface tensions were used in the calculations; we approxi- 
mated ■Jw-iiilJ-) ~ 7™2i(Mco) and 'yw2g{fJ') ~ 7™29(Mco). /J-pw 
denotes the chemical potential at pre- wetting, L = 00. In 
this approximation the thin-thick transition is independent 
of L. 



transition from a state with thin liquid films at both walls 
to the condensed state. There is a triple point where the 
prewetting transition intersects the condensation transi- 
tion line. The genesis of such a triple point was described 
in earlier pagers based on a simple DFT [2^ and a lattice 
gas model ,3C| . 

Now we return to the asymmetrically confined fluid. 
We have already seen that for T above the wetting and 
drying temperatures of the two walls, when L is large and 
the chemical potential is close to bulk coexistence, the slit 
contains a region of gas phase next to the drying wall and 
a region of liquid phase next to the wetting wall with 
a liquid-gas interface somewhere near the centre. We 
shall refer to this phase as the 'delocalised interface state' 
(i). Let us consider the prewetting transition that occurs 
at the solvophilic wall in the semi-infinite system. For 
dfi < Sfipu), apart from a thin liquid film at the solvophilic 
wall, the system contains mainly fluid at the bulk gas 
density. The situation in the finite L system at dfi < 
Sfipw is expected to be quite similar — the presence of the 
solvophobic wall will not alter the fluid density near to 
the solvophilic wall. The excess grand potential of this 
phase, which we shall refer to as the evaporated state (g), 
is given by 



when Sfi < and 



0{L- 



(43) 



when (5/i > 0, where ^W2g{^) is the surface tension of 
the thin film state in the semi- infinite fluid at wall2. 
(Note that when 5p, > Sfipw this wall2-thin film state 
is metastable). The bg/L^ terms in Eqs. (|42|) and 
arise from the interaction between the two wall- 
gas interfaces; the coefficient for this term is given by 
bg = {Pwi + Pw2 ~ /Og)Pg7re(7^/3. There is an extra 
term, 5p{pi — Pg)L, in the excess grand potential of the 
evaporated state (^5]) for Sp > because the bulk gas 
phase is metastable. At the prewetting transition in the 
semi-infinite system the wetting film thickness increases 
discontinuously so that the density profile has separate 
wall2-liquid and liquid-gas interfaces. In the asymmetric 
system the corresponding state has a thick film of liquid 
at wall2 separated from the gaseous region next to walli 
by a liquid-gas interface. This is what we term the de- 
localised interface state (i). The excess grand potential 
of this state (see Section III B 21) is what we calculated 
already 



^Ixi^eq', L, p) 
A 



Iwigip) + lw2l{fJ') + llg 
1 

'1?' 



;i:[p5p{L/af]+0{L-^), (44) 



where Yi[j55 p{L / aY] is a scaling function, given by Eqs. 
(1^^ and (|23p . Coexistence between the evaporated state 
and the delocalised interface state occurs when 



Pe 



(45) 



where Pevap is the chemical potential at evaporation. The 
phase diagram as a function oi 5p, the chemical potential 
deviation from bulk coexistence, and inverse wall sepa- 
ration L^^ is displayed in Fig. S) At large L evapora- 
tion occurs close to the chemical potential for prewet- 
ting at the isolated solvophilic wall, i.e., as L — > 00, 
Spevap — ^ Sppw and in Sec. IIVI we shall find the coex- 
isting density profiles near to the solvophilic wall at the 
evaporation transition are very similar to the thin and 
thick film density profiles for the semi-infinite fluid at 
the solvophilic wall. One observes in Fig. |4] that as L 
is decreased this transition moves to higher chemical po- 
tential, i.e., is closer to bulk coexistence, as the interac- 
tion of the liquid-gas interface with the solvophobic wall 
becomes more signiflcant. Evaporation occurs at bulk 
liquid-gas coexistence {6p = 0) for a wall separation of 



evap 



1/3 , ,1/3 



6,3 - bn 



lw2g(t^co) - lw2l(.fJ'Co) - llg 



1/2 



(46) 



^lx{L,p) 



Iwigip) + lw2g{y) + ^2 + '^('^ ) (42) liquid side of bulk coexistence — as for the evaporation 



For L < L'^yap the evaporation transition occurs on the 
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FIG. 4. Phase diagram for a fluid confined between asym- 
metric waUs at r = O.STc- The wall-fluid interparticle po- 
tentials are given by Eq. (|A2|) with pwi^wif = — l.lOSea"^ 
and Pw2^w2f ~ 1.723e(T~^. These parameters were chosen to 
ensure complete drying at isolated walli and complete wet- 
ting at walh at a temperature of T = O.STc- They ensure 
bi — &2 at this temperature. Eq. (|45|) was used to calcu- 
late the coexistence line between the two phases. The sharp- 
kink values were assumed for the coefficients bi, 62, 63 and 
bg. DFT was used to calculate the surface tensions for the 
liquid-gas jig, wall2-liquid 7™2i(/ico) and wall2-gas 'y^^gip.co) 
interfaces. Note that the walb-gas surface tension 7^29 (a'co) 
is the excess grand potential per unit area of the thin film 
state and refers to a metastable state at bulk coexistence fico- 
Again, we used the approximations "/w2i{fJ') ~ 7iu2i(Mco) and 
"fw2g{l^) ~ "fw2gil^co)- fJ-p-w denotes the chemical potential for 
prewetting at a single solvophilic wall2. 



transition in the slit with identical solvophobic walls (see 
Fig. [2]). However, for a given L evaporation occurs at a 
much lower value of (5/i in the asymmetric case. 

In the asymmetric system that we have studied using 
DFT (see Section HVl below) the solvophilic wall exhibits 
a first order wetting transition and, therefore, prewetting 
below the temperature of our investigations (T — O.STq) 
whereas the solvophobic wall is purely repulsive and is 
therefore dry at all temperatures. However, if the solvo- 
phobic wall potential was constructed such that there 
was a first-order drying transition at some temperature 
Td then a condensation transition, similar to the evap- 
oration transition described above and connected to the 
predrying transition in the semi-infinite fiuid at a sin- 
gle solvophobic wall, would be possible for temperatures 
To < T < (where T^^ is the predrying critical tem- 
perature). 



III. SUM RULES FOR THE CONFINED FLUID 



Recall that the total external potential acting on the 
fluid due to the two walls is 



and consider 
1 [dL 



dz 



sn 



^[^l-V{z;L)] 



6[^l-V{z■,L)]y^ 

^[^i-V{z:L)] 



dL 



dz pl{z) 



dL 



dzpLiz)VUL-z) 



(48) 



where we have used the following result for the functional 
derivative of the grand potential w.r.t. an external po- 
tential Vex{z)'- 



6n 



6[^i-V,,iz)] 



-P{z) 



(49) 



Pl{z) in (HS]) is the density profile of the fluid between the 
two walls at separation L. Let Vw2{z) = Vh{z) -\-Vatt2{z) 
where the first term is the potential due to a hard wall 
at z = and the second term is the attractive Lennard- 
Jones part of the wall potential. Then one can easily 
show 

ly'dL] =~I^~^P^'^-L^ I ^ pl{z)V^u2{L- z) 

(50) 

where Pw2,L is the contact density of the confined fluid 
at wall2 for wall separation L. Exchanging the two walls 
does not make any difference to the grand potential so 
we also have: 

^(^^j^ = -rV..i,L + J^^ dzpLiz)V:,,,iz). (51) 

At an isolated single hard wall we have the well-known 
result relating the pressure p of the fluid in the reservoir, 
far from the walls, to the density proflle: 



0+ 



dz Pociz)V^tt{z) 



(52) 



where Pw, oc is the contact density at the single wall and 
Poo (2) is the semi-inflnite density proflle at that wall. 
Substituting ([5^ and either (IHH) or ([51]) into the deflni- 
tion for the solvation force Eq. ([5^ we obtain: 



Pfs{L) — Pwl,L — Pwl,oo — /? 



Pw2,L — Ptu2,< 



dz pL{z)V^t^^{z) 
dz poo{z)V^tti{z) 



dz pl{z)V^u2{L ~ z) 



dz Poo{z)V^tt^{z) 



(53) 



(54) 



Viz;L) = V^i{z) + V^,2iL-z), 



(47) 



These results for the solvation force can be regarded as 
sum rules. In an exact treatment computing fs{L) from 
(1551) or must give identical results to those obtained 
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from evaluating the derivative of the excess grand poten- 
tial, i.e., dsn or ([22]). Non-local DFT treatments satisfy 
these sum rules. Thus within the context of DFT, the 
sum rules provide a stringent test of the accuracy of cal- 
culations. Provided we know the density profile of the 
fluid at one of the single walls then Eq. (|53p or Eq. (|54|) 
allows us to calculate the solvation force using the density 
profile for the confined fluid at just one value of the wall 
separation L rather than performing the derivative in 
Eq. p2p numerically. Comparing results from both Eqs. 
([5^ and ([51]) provides a consistency check for our DFT 
calculations below fSec. IIV|) . Such checks are important 
given the small free energy differences and therefore the 
resulting sensitivity of the phase transitions and scaling 
functions to the numerics. 



IV. DFT RESULTS 

In this section we present numerical results, obtained 
using a fully microscopic DFT approach, for a fluid con- 
flned between planar opposing walls. The fluid-fluid and 
wall-fluid interaction potentials are described earlier in 
Section III Al The excess hard sphere part of the free en- 
ergy functional J-"^^ was treated by means of Rosenfeld's 
fundamental measures theory (3T| and the attractive part 
of the fluid-fluid interaction potential was treated in 
mean-field fashion. The functional used was the same 
as in Refs. (32l . [33| where it was used to describe drying 
at the surface of a sphere. The grand potential functional 



dridr2 p(ri)p(r2)0 att U 1 



r2|) 



+ / p{v){V{v)^^J)Av. 



(55) 



where the density profile p(r) = p[z) and the external 
potential is given by (gT]): V{t) = V{z) = V{z]L). Fid[p\ 
is the Helmholtz free energy functional for the ideal gas. 
4>att is given in ([1]). The equilibrium density profile was 
found by minimising the grand potential functional 
and the corresponding equilibrium excess grand potential 
r^ea; was calculated. The numerical methods employed to 
obtain the equilibrium density profile are described in 
detail in Chapter 3 of Ref. 34 1 . 

The bulk coexistence curve calculated from the func- 
tional ([55)1 has the standard mean-field form with the 
critical temperature /csTc/e = 1.415 and the critical den- 
sity pccr^ = 0.2457. All our results were obtained at the 
same temperature T — O.STp. Our choice of parame- 
ters, Pwi^wif = — l.lOSeo-^"^ and Pw2^w2f = 1.723ea'~'^, 
ensures complete drying at an isolated walli and com- 
plete wetting at wall2 and also results in equal coeflScients 
bi — 62 in the effective interfacial potential (see Eqs. [U 
[S]and[7]) at this temperature. 

Figure [5] shows density profiles of the confined fluid, 
with wall separation L = lOOti, at liquid-gas coexis- 
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FIG. 5. Density profiles p{z) at various chemical potentials 
for the asymmetric slit with a wall separation of L = IOOct 
and temperature T — O.STc- The full line shows the profile at 
bulk liquid-gas coexistence, Sfi — 0. Because bi = 62 the gas- 
liquid interface is located at the mid-point between the walls 
for Sp = 0. Profiles a-d correspond to chemical potentials 
on the liquid side of coexistence, f^Sp, = 10~®, 10~^, 10"* 
and 10""^ respectively. Profiles e-h occur on the gas side of 



coexistence for /35/i 



-10" 



-lO^" and -10" 



fence and also for the fluid on both sides of coexistence. 
The walls are chosen to be 'antisymmetric' in that the 
leading order term in the interaction potentials between 
each wall-fluid interface and the liquid-gas interface is the 
same (i.e., 61 = 62). Consequently the liquid-gas interface 
lies in the centre of the slit at bulk coexistence 5p — 0. 
If we increase the reservoir chemical potential away from 
coexistence into the liquid region of the phase diagram 
so that 6p is positive then the gas-liquid interface moves 
away from the centre and towards the drying wall, de- 
creasing the volume of the now metastable gas phase. If 
Sp is negative so that the reservoir is in the stable gas 
phase then the volume of liquid phase shrinks and the 
gas-liquid interface is closer to the wetting wall. The 
shape of the density proflles near to the walls is largely 
unaffected by these changes in the chemical potential. 
Packing effects in the liquid adsorbed at the wetting wall 
result in oscillations in the density proflle in the close 
vicinity of the wall. The density profile is monotonically 
decreasing on approaching the drying wall, with a con- 
tact density slightly below the bulk gas value [32l|. 

In order to test the effective potential prediction for 
the film thickness scaling function, A [l36p{L/(T)^'\ , Eqs. 
(fT5| and (fTH]) . we calculated the excess adsorption per 
unit area ([T7]) at various chemical potentials and for four 
different sized systems, L/a = 500, 250, 100 and 50. The 
thickness of the gas film l^q is related to the excess ad- 
sorption r by Eq. (fT8|) or Eq. (IT9|) depending on whether 
the chemical potential of the reservoir is on the liquid or 
the gas side of bulk coexistence. The coefficients 61 and 
62 appearing in the excess grand potential for the system 
(Uni and adsorption scaling function ([16]) are independent 
of the precise definition for Zgg [11] but the next to lead- 
ing order terms with coefficients ci and C2 are affected 
by this choice. Figure [H] displays the DFT results and 
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FIG. 6. Plot to illustrate adsorption scaling in the asymmet- 
ric slit on the liquid side of bulk coexistence, i.e., Sn > 0, 
where the adsorption is related to the thickness of the gas 
film by leq = —r/[A{pi — Pg)]. The full line shows the pre- 
diction of the effective potential approach, i.e., the function 
-A [l35ii{L/af] valid as L ^ oo, given by Eq. (|16|l . The 
symbols are DFT results for various wall separations. 



the effective potential prediction, Eqs. dTSl) and p^ . for 
tfie adsorption scaling function on the liquid side of bulk 
coexistence. There is good agreement between the DFT 
results and the scaling prediction over a very large range 
of values of {L/a)^f35fi and we found equally good agree- 
ment on the gas side of bulk coexistence (not shown). 
Near to bulk coexistence, i.e., for small values of \6ii\, 
the liquid-gas interface is close to the mid-point between 
the two walls so that |r|/[A(/9; — Pg)L] « 0.5. The inter- 
face begins to move away from the centre towards one of 
the walls as \5p\ is increased and the term 5p{pi — pg)l 
in the grand potential ((T5)) becomes significant at around 
{L / aY (55^ = 10. Eventually, for large chemical potential 
difference the liquid-gas interface is close to one of 
the walls and the terms in either or [L — l)^^ in the 
grand potential (jl3|) become important and the approx- 
imation which led to our scaling prediction (jl5l and I16p 
breaks down. This can be gleaned from our DFT results 
for L/cr = 50 and L/a = 100 which deviate slightly from 
the effective potential prediction at some of the larger 
values of {L/aYl3\5p\. 

In Section III B 21 we showed that the surface tension 
of the system can be expressed as the sum of the surface 
tensions of the three individual interfaces (walli-gas, gas- 
liquid and wall2-liquid) plus a term proportional to 
multiplied by a scaling function Yi[j35ii{L / a)^], see Eq. 
(I^TI) . In Figure [7] we compare the predicted scaling func- 
tion (1^^ with DFT results obtained on the liquid side of 
bulk coexistence. For this case and on the gas side of bulk 
coexistence (not shown) there is good agreement. The 
very small discrepancy between the DFT results and the 
prediction from Eq. ([^^ at small values of {L/aYP\5p\ 
probably arises because the difference between the sur- 
face tensions, 7(i,/z) - jwigiPco) - lw2i{p co) llgi IS 
very small compared to the values of the surface tensions 
themselves. When L is large any error is magnified when 



FIG. 7. Plot to illustrate surface tension scaling in the asym- 
metric slit on the liquid side of bulk coexistence 5p > 0. Sym- 
bols denote DFT results for various wall separations. The in- 
dependently calculated walli-gas and wall2-liquid surface ten- 
sions at coexistence, and the gas-liquid surface tension have 
been subtracted from the total surface tension of the system. 
When multiplied by and plotted against SpL^ the data 
for various different L and 5/i collapse onto a single curve. 
The full line is the scaling function /?E \j35p,{L / a)"^^ , valid for 
L ^ oo, calculated from Eqs. (pT|) and ((22)) . 



multiplied by , which explains why the greatest devia- 
tion from the theory is for the DFT results at L = 500cr. 

There are three different ways that we can extract the 
solvation force fs{L) from our numerical data. It can be 
calculated, using Eq. ([55]). from the change in surface ten- 
sion S'j when the distance between the walls is increased 
by a small increment SL, i.e., fs{L) = —(6j/SL)fj_^T,A- 
Alternatively it can be found from the density profile by 
applying a sum rule at either of the walls, as described 
in Section IIIII In Figure [8] we plot the solvation force, 
obtained via these three methods using wall separations 
of L = IOOct and L — lOlcr, as a function of the chem- 
ical potential deviation from bulk coexistence. There is 
very good agreement between the values for the solvation 
force obtained using sum rules at the two different walls 
for each value of L. The values calculated from the dif- 
ference in the surface tension between the systems with 
L — 100(7 and those with L = lOltr lie in between the 
sum rule results from the two separate systems. Having 
confirmed that the three numerical methods for obtaining 
fs are consistent we chose the most convenient to test the 
scaling behaviour predicted by Eqs. and ([55]) . For 
systems on the liquid side of coexistence the solvation 
force was obtained most easily using the sum rule at the 
wet wall ([54]). DFT results for at L = 500cr are close 
to the predicted scaling function but there is more dis- 
crepancy as L is decreased and for L — 50(T there is a 
significant deviation — see Fig. (S) This can be explained 
by the presence of higher order terms (proportional to 
L~^) in the excess grand potential which we neglected in 
our derivation of Eqs. ([M)) and ([55)) . If these terms are 
included then it is no longer possible to express the solva- 
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FIG. 8. DFT results for the solvation force in the asymmetric 
slit as a function of the deviation from coexistence Sjj,, for 
5fi > 0. The circles (•) are data obtained by directly measur- 
ing the change in the surface tension between L = lOOcr and 
L = 101a. The other symbols are results acquired using the 
contact density sum rules at the two walls for L — lOOcr and 
L = lOla— see Eqs. ^ and ((54]). 
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FIG. 9. DFT results for the solvation force in the asymmetric 
slit multiplied by L^, obtained using the sum rule (|54|) at the 
wet wall (wall2) for Sjj, > 0. The full line is the predicted 
scaling function f^Fg [/35/i(L/(T)'^] , valid as L — > oo (|35|l . The 
dashed line is the prediction for L — lOOcr which includes 
sharp-kink terms of order in the effective potential. 



tion force in terms of a scaling function in the manner of 
Eq. In Fig.iniwe have plotted the effective potential 
prediction for the solvation force for L = lOOcr including 
the next to leading order contributions with their sharp- 
kink coefficients. The agreement with the DFT results 
for L = 100(7 is surprisingly good as we expected the 
coefficients ci, C2 and C3, Eqs. ©, pU)) and ([TT|) to be 
affected by the details of the density profile. 

The effective potential approach predicts a suscepti- 
bility which diverges as at bulk liquid-gas coexistence 
(PS)) . Results for the total susceptibility, x, are shown 
in Figure [To] for 5fi > 0. The magnitude of the suscep- 
tibility is greatest at bulk coexistence Sfj, = when the 
liquid-gas interface lies at the mid-point between the two 
walls; in this position the external potential exerted on 



FIG. 10. Plot to illustrate susceptibility scaling in the asym- 
metric slit for > 0. The symbols give DFT results for 
various L and Sjj,, obtained by finding the difference in the 
adsorption AT between results at chemical potentials differ- 
ing by f^Afi = 1 X 10~®. The full line gives the predicted 
scaling function /3"V®C [/3(5/i(L/cr)^] from Eq. (gZ)), valid as 
L — )■ 00 and the dashed line is a prediction which includes 
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the fluid by the two walls is very weak and the free energy 
cost incurred by shifting the interface is smallest. As \Sjj,\ 
increases the magnitude of x begins to fall sharply as the 
liquid-gas interface moves away from the centre of the 
slit— at around {L/aY/SSfj. = 10. The DFT resuhs are 
in agreement with the predicted scaling function (|27p for 
L — 500(7 but begin to deviate for smaller L. Including 
higher order terms (those proportional to {L — 
and in the effective potential (1X31) . with sharp-kink 
coefficients) does account for some of the disparity, as is 
illustrated by the prediction for L = lOOcr (dashed line 
in Fig. 1101) which is much closer to the corresponding 
DFT results. The remaining discrepancy is probably due 
to the coefficients of these terms being modified because 
the DFT density profiles are smooth rather than sharp- 
kink-like. 

Figure [TT] shows the local susceptibility x(2;-^, de- 
fined by Eq. ([50]) . i.e., the change in the density pro- 
file p{z) with chemical potential, for a wall separation 
of 250ct, evaluated at bulk liquid-gas coexistence fico- 
The peak in the magnitude of the susceptibility coin- 
cides with the position of the liquid-gas interface lying 
at the centre of the slit confirming that the most sig- 
nificant change in the density profile as the chemical 
potential is varied involves the displacement of this in- 
terface. This conclusion is backed up by the similar- 
ity between the local susceptibility —{dp{z)/dfi)L,T and 
the product {dleq/dp)L.T{dp{z)/dz), see Eq. ([5T|) . which 
is also plotted in Fig. [TTJ (The quantity (dleq/ dji) j, 
was calculated from the change in the adsorption using 

{di^jdp)^., = -[dT/dtj)L,T[A{pi - Pg)]-^-) 

In Section lll Cl we predicted that a first order capillary 
evaporation transition could occur in our asymmetric slit 
and used the effective potential approach to calculate the 
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FIG. 11. The magnitude of the local susceptibility —xi^) = 
— X(«;i, A*) = idp{z)/dfi)j^ J, near the midpoint between the 
two walls obtained by subtracting profiles at — Pjico and 
= jl^co + 1 X 10~^" for a wall separation of L = 250ct (■). 
The line gives the product ~{dleq/dfi)L,T{dp{z) /dz). 
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FIG. 12. The coexisting density profiles of the 'delocalised 
interface state' (full line) and the evaporated state (dashed 
line) for I/eSap = 11.4(t and Sjj, = 0. 



chemical potential at this transition as a function of wall 
separation (Fig. |4]). The line of evaporation transitions 
crosses bulk liquid-gas coexistence at a predicted value 
of L'^yap = 13cr, obtained from Eq. (j^Hl) using sharp-kink 
values for the coefficients 6i, 62, 63 and bg and DFT re- 
sults for the surface tensions Jw2gil^co), Iw^iifJ^co) and 
"fig. (Note once again ^^^gifJ-co), is the excess grand po- 
tential of the wall2-thin film state which is metastable at 
A* = Mco-) Comparing the excess grand potential (calcu- 
lated from DFT density profiles in the asymmetric slit) 
for the evaporated state with that of the state with that 
of the 'delocalized interface state', as a function of L gave 
the slightly smaller value of Ll°ap = 11.4a. It is not sur- 
prising that the effective potential prediction of L'^°ap is 
not particularly accurate at such small wall separations. 
Figure [T^ shows the density profiles of the two coexisting 
states at L'^°ap and it is clear that the density profile for 
the delocalised interface state is far from the sharp-kink 
approximation for the profile which forms the basis of the 
effective potential treatment in Section |TT1 

At larger wall separations, L > L'^°apy capillary evap- 
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FIG. 13. The coexisting density profiles of the 'delocalised in- 
terface state' (full line) and the evaporated state (dashed line) 
for L = lOOcr and pSfievap = -4.1344 x 10"^. The profiles 
are almost identical to those corresponding to the coexisting 
thick and thin films at the prewetting transition for a single 
solvophilic wall2. 

oration occurs when the system is on the gas side of 
bulk coexistence, i.e., Sfievap < (see Fig. |4]). This is 
in contrast to the situation for a system with identical 
drying walls where capillary evaporation may occur only 
as bulk coexistence is approached from the liquid side 
(see Fig. [2]). For wall separations L > 50 evaporation 
is predicted to occur at chemical potentials very close to 
SfJ'pw, the chemical potential at the prewetting transition 
at the isolated solvophilic wall. Using DFT we found that 
PSfipw — —4.1368 X 10~^ for prewetting at (single) wall2, 
c.f., the value —3.4 x 10"'^ from the effective potential 
approach used in Fig. [51 The chemical potential at the 
evaporation transition for a wall separation oi L — lOOcr 
was found to be PSfXevap = —4.1344 x 10^"^ which is in- 
deed very close to S^pw Figure [T3] displays the coexisting 
density profiles at capillary evaporation for L = lOOcr. 
Unlike the profiles shown in Fig. [S] for smaller values of 
|(5/i|, at this chemical potential the gas-liquid interface 
in the 'delocalised interface state' has moved away from 
the centre of the slit and closer to the solvophilic wall2. 
The evaporated state is filled with gas apart from a thin 
film (thickness sa la) of higher density fiuid next to the 
solvophilic wall. The two density profiles in Fig. [T^ are, 
apart from the region very close to the solvophobic wall, 
almost identical to those of the coexisting thick and thin 
films at the prewetting transition at wall2. 

V. DISCUSSION 

In this paper we made a detailed investigation of the 
'delocalised interface phase' or 'soft mode phase', identi- 
fied originally by Parry and Evans @ for a model fluid de- 
scribed by long-range (— r~^) dispersion potentials which 
is confined between competing planar substrates. We 
summarize our main results as follows: 

1. Using an effective potential approximation (Sec- 
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tion In]) we derived scaling forms for the Gibbs 
adsorption T{L,^), the surface tension 7(L,/i), 
the total susceptibility x{L, m) and the solvation 
force fs that depend on the dimensionless prod- 
uct {L/a)^/3Sfi, where the power of three reflects 
directly the power of —6 in the intcrparticlc poten- 
tials. The scaling forms were confirmed by numer- 
ical results from calculations based on a non-local 
DFT (Figs.il[ll[9]and[l0l). Provided the waU sep- 
aration L is sufficiently large, typically L > lOOa, 
there is good agreement with the explicit scaling 
predictions. Perhaps most striking is the predic- 
tion that x(L,/^co) i"* — see Eq. ([29|. This re- 
sult corresponds to a transverse correlation length 
that diverges as ^|| ~ L^, for L — > oo Both 
these results imply that the central liquid-gas in- 
terface undergoes extremely large interfacial fiuc- 
tuations, correlated over very long distances paral- 
lel to the walls. These interfacial fiuctuations re- 
flect the smallness of the free-energy cost of shift- 
ing the whole liquid-gas interface and are charac- 
teristic of the 'soft mode phase'. Recall that for 
S/j, = the liquid-gas interface is constrained by 
only the tails of the interfacial potentials which de- 
cay as a power-law for dispersion forces). By 
contrast, for short range forces the decay is expo- 
nential which leads to the result quoted in Section 
H x{L,fico) - e«^/2, for L ^ oo 9]. 

2. The central liquid-gas interface is subject to an ef- 
fective repulsive force from each wall as walli en- 
deavours to become dry and wall2 to wet. The re- 
sulting solvation force is repulsive. For dispersion 
forces we find fs ~ +L~'^ for Sfj, = 0. Specifically, 
in the perfectly antisymmetric case with bi = 62, 
the liquid-gas interface lies in the centre of the slit 
and from Eq. (I35|) 

/,(i,Aico) = -^(1652 + 263). (56) 

This result should be contrasted with that for 
short-range forces in the 'soft mode phase' where 
fs ~ -he-''^/^ for L — ^ c» 0. Employing the 
sum rules of Section IIIII we calculated the solva- 
tion forces from DFT using three different routes, 
demonstrating their consistency and the accuracy 
of our numerics — see Fig. [S] 

It is interesting to compare our results for fs with 
those obtained using an equivalent DFT approach 
for a different model fluid conflned between two 
identical walls [S^ . In this study the wall-fluid po- 
tentials decayed as and the fluid-fluid pair po- 
tentials were short-ranged, i.e., truncated Lennard- 
Jones. Away from any phase transitions and bulk 
criticality the solvation force was found to decay 
as fs ^ +L~^- However, including the full —r~^ 
tails of the fluid-fluid pair potential results in an 
additional attractive —L~^ contribution to fs [36| . 



The asymptotic decay of the solvation force for a 
system with identical walls can be either attractive 
or repulsive depending on the choice of parameters 
in the (power-law) wall-fluid and fluid-fluid poten- 
tials and on the thermodynamic state point. 

3. In Section III CI we used the effective interfacial po- 
tential approximation to investigate possible phase 
behaviour for our model fluid confined between 
asymmetric walls. Figure |4| shows the predicted 
phase diagram for the single temperature T = O.STc 
that we studied using DFT in Section IIVI Depend- 
ing on the wall separation L capillary evaporation, 
i.e., a transition to a dilute gas state with only a 
thin adsorbed film of liquid-like density at wall2, 
can occur not only for (5/i > (the situation that 
pertains for capillary evaporation in the case of 
identical walls) but also for 6fi < 0. Fig. IT2| shows 
the coexisting density profiles calculated from DFT 
when capillary evaporation occurs at Sfi = 0, i.e., 
for the wall separation L'^yap = 11.4cr. Evapora- 
tion occurs because for smaller wall separations, 
L < L'^°apj the free energy cost from the inter- 
face interactions plus the surface tension of the cen- 
tral liquid-gas interface is so great that it is more 
favourable for the slit to be filled by the evapo- 
rated 'gas', despite the fact that the surface ten- 
sion Jw2g(p) between gas and (solvophilic) wall2 is 
large. For L > L'^°ap capillary evaporation occurs 
on the gas side of bulk coexistence, i.e., Sfj. < 0. 
However, in this regime for large L the evapora- 
tion we observe is closely connected to the prewet- 
ting transition that occurs at an isolated {L = 00) 
solvophilic wall. Thus for L = lOOcr DFT yields 
l3{Sfievap — SfJ-pw) = 2 X 10^^ and Fig. IT^ displays 
the coexisting density profiles at this evaporation 
transition. These profiles are almost identical to 
those at the prewetting transition. Our present 
study describes evaporation, rather than condensa- 
tion, because the solvophilic wall2 undergoes a first 
order wetting transition at Tw < 0.8Tc whereas 
the purely repulsive solvophobic walli is dry at all 
temperatures between the bulk triple and critical 
temperatures; there is no pre-drying transition. In 
a perfectly antisymmetric system, such as an Ising 
model subject to equal but opposite surface fields 
(hu = —hi) for which walli undergoes a first or- 
der 'drying' transition at the same temperature as 
wall/) undergoes the equivalent wetting transition, 
the evaporated state with a thin film of 'liquid' at 
wall/) would coexist with the condensed state with 
a thin film of gas at walli and with the 'delocalised 
interface state' at bulk coexistence Sfi = 0, result- 
ing in a triple point at L = = The 
admirable review by Binder et.al. [3] provides an 
illuminating summary of the genesis of such triple 
points in Ising models, focusing mainly on temper- 
ature dependence at fixed L rather than fixing T 
and varying L as we do here. 
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Miiller et.al (see Ref. p?] for a summary) also found 
triple points in their self-consistent field treatments 
of symmetric AB binary polymer blends confined 
between perfectly antisymmetric walls, i.e., one 
wall attracts the A component with exactly the 
same strength as the other wall attracts compo- 
nent B and there is a first-order wetting transition. 
Although the interparticle potentials in these stud- 
ies are short-ranged it is probable than many of 
the gross features of the phase diagram found by 
Miiller et.al. will also apply to our system. 

Our results pertain to non-retarded (London) dispersion 
forces: we deliberately chose to consider interparticle po- 
tentials that decay as —r~^ since these describe the mod- 
els that are conventionally used in the physics of liquids. 
In the physics of wetting it is well-known that phase tran- 
sitions at a single wall (L = oo) depend profoundly on 
whether the interparticle pot entials are short-ranged or 
long-ranged (power-law) 27]. Indeed the critical expo- 
nents describing continuous wetting transitions depend 
on the specific power. Relevant to the present study is 
the observation that for fluids subject to wall-fluid po- 
tentials decaying as —z~^"'~^^\ corresponding to a wall 
particle-fluid particle potential decaying as — r~("+'*\ the 
leading term in the excess grand potential ^ goes as 
where n = 2 for non-retarded dispersion forces. For the 
complete wetting transition (from off bulk coexistence, 
with (5/i — 0) the upper critical dimension for systems 
with power-law potentials is given by [1, [l^l : 

d* = 3-4(n + 2)-i (57) 

Thus for complete wetting in the presence of non- 
retarded dispersion forces d* = 2 and mean-field pre- 
dictions for the critical exponents will remain valid for 
all dimensions d > 2. We recall now the heuristic scaling 
ansatz in Section 4 of Parry and Evans Q • This argues 
that the only relevant scaling argument in the 'soft mode 
phase' must be the ratio L/{2leq), where the equilibrium 
thickness of the (drying) film is given by the complete 
wetting resuh leg ~ I^A'I"'^" with = {n + 1)"^ for 
power-law potentials. The idea is that the maximum 
thickness of the film at equilibrium is L/2, at 6fj, = 0, 
and that the film thins according to the complete wet- 
ting behaviour. Note that there is no temperature scaling 
variable: it is assumed that the temperature is sufficiently 
far below Tc and sufficiently far above the interface de- 
localization temperature that it is irrelevant. Only the 
chemical potential deviation Sfi is relevant for scaling. 
The resulting scaling predictions ^ for a general power- 
law potential are x{L,Hco) and ^|| ^ 
for dfj, = and L — ;> oo.The solvation force is predicted 
to decay as -l-L"*^""'"-'^) in this 'soft-mode' regime. 

Our present results, based on the effective interfacial 
potential treatment of Section [II] and the DFT calcu- 
lations of Section IIVI are evidently consistent with the 
scaling ansatz for the particular case n — 2. Clearly both 
are mean-field treatments; they both fail to capture all 



of the effects of capillary- wave fluctuations. Nevertheless, 
if the heuristic scaling ansatz of Ref. @ is correct, and 
there seems no reason to doubt its validity, and given 
Eq. (1571) for the upper critical dimension for complete 
wetting, it follows that for all power-law potentials (n fi- 
nite) d* < 3 so that the mean-field predictions should be 
valid as regards the dependence of thermodynamic quan- 
tities and the correlation length on the wall separation 
L, for L — > oo. This observation is important since it 
implies that the leading order scaling functions that we 
derived from the (mean-field) effective interfacial poten- 
tial with n = 2, and checked using the microscopic DFT 
approach, should remain valid in the presence of fluctu- 
ations. Moreover, so should the corresponding analysis 
for any finite n. Of course, in an exact treatment the 
coexisting bulk densities pi and pg entering the Hamaker 
constants bi, &2 and 63 — see Eqs. ([6l[7]and[8]) — should be 
replaced by their exact counterparts for the model fluid 
under consideration. 

What is omitted in our mean-field DFT approach? 
Capillary-wave fluctuations must broaden the central 
liquid-gas interfacial density profile in the 'soft mode 
phase'. This effect is not captured within the simple DFT 
employed here, i.e., the results shown in Fig. [5] showing 
rather sharp interfaces do not encompass such broaden- 
ing. Computer simulations for the same model would 
exhibit broader interfacial density profiles. In spatial di- 
mension d = 3 the extent of the capillary wave broaden- 
ing is given by S,±_ ^ (ln^||)^/^. For the large wall separa- 
tions that we consider this broadening is not insubstan- 
tial and is certainly on the scale of the intrinsic width 
that emerges from our DFT calculations. Once again we 
contrast this result with the case of short-ranged forces 
where ^|| - e''-^/'* [9] and thus - L^/^, for Sp = and 
L — > 00. This last result was confirmed (for relatively 
small L) by measurements of the width of the magnetiza- 
tion profile in Ising model simulations, e.g., the summary 
in Section 3.4 of Ref. [|. 

We emphasize that the analysis in the present paper is 
for a single temperature that lies in the 'delocalised in- 
terface phase'. Unlike earlier studies that focused on sys- 
tems with short-ranged forces 13-15i| we have not 
attempted to investigate the localisation-delocalisation 
transition that is intimately linked to the wetting transi- 
tion at an (isolated) confining wall. In the original stud- 
ies Q the wetting transition was continuous (critical) 
and there were detailed predictions for the location of 
the critical point and the nature of the criticality in the 
confined system with perfectly antisymmetric walls. Sub- 
sequently these were confirmed by results from Ising sim- 
ulations d, [isl - fTsI . For perfectly antisymmetric walls the 
critical temperature Tql (of the interface delocalization- 
localisation transition) lies below the single wall critical 
wetting temperature Tw for any finite L and in the limit 
L — >■ 00, TcL — >■ Tw- In the present analysis we de- 
liberately chose the wall-fluid interaction parameters in 
order that the walls were antisymmetric to leading or- 
der in the excess grand potential ([2]). Speciflcally we set 
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6i, the Hamaker constant for drying at walli to be equal 
to 62, the Hamaker constant for wetting at wall2. Of 
course, the Hamaker constants depend on temperature 
through their dependence on the bulk coexisting densi- 
ties, see Eqs. (|6]) and ([7]). 61 and 62 are equal at only one 
temperature. At a lower temperature these quantities are 
not identical and the system is not antisymmetric. This 
makes locating the localization-delocalisation transition 
rather demanding. 

We have also deliberately avoided the bulk critical 
regime, near Tc, where the bulk correlation length S,b ^ 
L. For competing walls we expect a universal 'critical 
Casimir' solvation force that decays as fs ~ +L~'^ ex- 
actly at criticality. Thus for d=3 the power-law decay of 
the critical Casimir force is the same as that found for 
non-retarded dispersion forces in the 'soft mode phase' 
below Tc- The amplitudes are, of course, different. Ascer- 
taining the precise nature of cross-over as T approaches 
Tc from below is challenging. 

We should comment on studies of other model fluids in 
asymmetrical confinement. As we mentioned in Section 
m colloid-polymer mixtures are a particularly attractive 
system to study both from an experimental and a the- 
ory/simulation perspective — especially if one adopts the 
simplest Asakura-Oosawa-Vrij (AO) model for these bi- 
nary mixtures Q . Within the context of this fluid model, 
a hard-wall will favour wetting by the colloid-rich phase, 
as a result of depletion interactions, while a suitably tai- 
lored wall, with appropriate coating, can be made to 
favour the polymer rich phase @. Since depletion inter- 
actions arise from excluded volume considerations they 
are short-ranged. Indeed they are strictly finite ranged 
in the case of the AO model. (Experimentally it is as- 
sumed there are no residual dispersion forces, i.e., refrac- 
tive index matching is perfect.) Nevertheless, several of 
the results obtained by De Virgiliis et.al. [l3j[3 in sim- 
ulations of the AO model are close to those that we find 
for our atomic model. For example, the density profiles 
for the colloid species in the 'delocalised interface phase' 
(see Figs. 2b and 15d, where L — 10 colloid diameters, of 
Ref. [iSl p are very similar to what we find in the corre- 
sponding 'delocalised interface state' (see Fig. [T^] of the 
present). It is expected that several features the phase 
behaviour ascertained for the AO model in asymmetric 
confinement will pertain to the present fluid model. Of 
course, the details of the decay (or divergence) of ther- 
modynamic functions will be different. 

Returning to the present one-component model one 
might argue that the choice of a purely repulsive sub- 
strate (walli) is inappropriate for any real fluid. In reality 
there is always some residual attractive dispersion inter- 
action between the fluid and the substrate and Young's 
contact angle is always < tt; one does not have com- 
plete drying. However, one can easily circumvent this 
objection by considering a binary mixture with dispersion 
forces. Similar to the case of polymer blends, one can en- 
visage a binary molecular mixture confined by asymmet- 
ric walls that favour one species or the other so that in a 



certain temperature range one wall can be completely wet 
by the phase rich in species A and the other completely 
wet by the phase rich in species B. It is straightforward 
to extend the present effective interfacial potential ap- 
proach to the binary case and evaluate the appropriate 
Hamaker constants that determine the scaling functions. 
The phenomenology is somewhat richer because we deal 
with a mixture but the basic scaling predictions in the 
'soft mode phase' remain the same. 

It is important to review the length scales that we con- 
sidered in this study. In order to confirm the scaling pre- 
dictions it was necessary to perform DFT calculations 
for extremely large wall separations — in some cases up 
to 250(7. Clearly this range is far beyond the current 
realms of molecular simulation. Since the 1930s it has 
been recognized that dispersion forces play a key role 
in the physics of wetting and in many other aspects of 
confinement. It is clear that a full understanding of sev- 
eral subtle aspects of wetting phenomena requires a mi- 
croscopic approach that incorporates these long-ranged 
forces and that treats short-ranged correlations, arising 
from packing of the 'particles', in a realistic fashion. Our 
present study and our earlier one on wetting/drying at 
a curved substrate [32] demonstrate that classical DFT 
provides a successful approach, albeit one that treats the 
attractive interactions at mean-field level. One should 
also note that retardation becomes relevant in real flu- 
ids and for real substrates. For sufficiently large inter- 
nuclear separation r the atom-atom pair potential crosses 
over from —r~^ to — decay. This implies cross over to 
ultimate — decay of the wall-fluid potential and to a 
leading contribution of in the excess grand potential 
([5]). Were we to attempt a more realistic treatment of 
interactions for large separations, where retarded forces 
are relevant, the ultimate power laws would be different 
from those obtained for non-retarded dispersion forces. 
For example, the solvation force decays ultimately as 
fs ~ +L~'^, rather than as +L^^. This observation is 
important since we showed in our microscopic calcula- 
tions for the non-retarded case (see Fig. ^ that agree- 
ment with the scaling limit is achieved accurately only 
for very large L, say > lOOcr. In real systems retardation 
will kick in at similar separations. 

Can the results of our present study be related to ex- 
periment? By measuring the force between the tip and 
a substrate in an atomic force microscope or between 
the crossed cylinders in a surface force apparatus (SFA) 
one can obtain information about the solvation force. In 
particular, by using the Derjaquin approximation, one 
can relate the force measured in the SFA to fs{L, ji), the 
quantity we calculated and which pertains to two (infi- 
nite area) parallel walls [37| . Our present study was mo- 
tivated in part by an experimental study 38] that used 
a SFA to investigate both the normal force and the re- 
sponse to shear deformation of water, in contact with a 
reservoir at normal temperature and pressure, confined 
between hydrophobic and hydrophilic surfaces. The ti- 
tle of Ref. 38] refers to a Janus Interface. In the ex- 
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periments the hydrophobic surface was either a cyhnder 
of mica coated with a self-assembled monolayer of oc- 
tadecyltriethoxysilane (OTE) or mica with a thin film of 
silver then coated with a self-assembled octadecanethiol 
layer. The hydrophilic surface was untreated mica which 
the authors assumed to be wet completely by water. It 
was observed that the shear response was extraordinarily 
noisy and indicated a distribution of relaxation processes. 
The authors mention 'giant fluctuations (of the dynam- 
ical shear re&ponses)' and allude to the work of Parry 
and Evans 0. Without being explicit they appear to 
infer that their experimental configuration corresponds 
to that of the 'soft-mode phase' in Ref. @, |1] and they 
write about a 'flickering, fluctuating complex' in which 
the capillary wave fluctuations are somehow thwarted. 

Two observations are relevant in the light of our 
present analysis of the soft-mode phase: 

(i) The measured static force-distance profile (see in- 
set to Fig. 1 in Ss]) is attractive over a wide 
range of surface separations, from 1000 down to 
about 10 molecular (water) diameters. This con- 
trasts sharply with what we calculate for fluid con- 
finement between competing solvophobic (complete 
drying) and solvophilic (complete wetting) walls. 
For this case we find the solvation force is repulsive 
at large separations throughout the region of the 
soft-mode phase. At small separations L, where 
the confined fluid is a dense liquid, the solvation 
force can oscillate as a function of L due to pack- 
ing effects. 

(ii) The hydrophobic surfaces prepared in [sll are, of 
course, not completely dry. The contact angle mea- 
sured for the OTE surface is 110 ± 2° and that for 
the thiol surface is 120 ± 2°. Thus the situation 
realized in the experiment js^ does not match that 
described in d, and which is considered in detail 
in this paper. It difficult to see why the particu- 
lar choice of wet and partially dry surfaces stud- 
ied in the SFA experiments [3^ could give rise to 
a wildly fluctuating liquid-gas interface — see also 
the comments of Pertsin and Grunze [s^. Nev- 
ertheless the experimental observations remain in- 
triguing and worthy of further investigation. 
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Appendix A: The wall— fluid potential 

We model the substrate as a block of (wall) particles 
of constant, uniform density pw The wall cannot be 



penetrated by the fluid particles, i.e., the wall-particle 
interaction potential is infinitely repulsive at the plane 
of contact. This hard-core repulsion between the parti- 
cles results in an excluded volume, of width dw, where 
the density vanishes, at the boundary between the wall 
and the fiuid, as illustrated in Fig. [TJ The wall is struc- 
tureless, i.e., parallel to the interface the density is taken 
to be uniform. The attractive potential exerted by the 
wall on the fluid particles can be calculated by first con- 
sidering the interaction potential between individual wall 
particles with fluid particles. This is chosen to be of the 
same form as the fluid-fluid inter-particle potential. The 
attractive part is given by 



Ct{(''l2) 








(?s) - 









^12 '^min,w f 
^12 '^min,w f 

(Al) 

(c.f. Eq. [1]) where ri2 is the distance between the wall 
particle and the fluid particle and e^,/ and a^jf are the 
strength and range of the wall-fluid potential, respec- 
tively and rjnm,wf = 2^/^(Tu,/. The total potential, ex- 
erted by the entire wall, on a single fluid particle at a 
perpendicular distance z from the contact plane (see Fig. 
[T]) is found by integrating the inter-particle potential over 
the semi-infinite volume of the wall: 



z <0, 



C+,n, d^' r 2nr^Z{iV^^TJ^) dr, z > 0. 



We obtain 



[PwVwf{z + dw), z > 0, 

where the function Vwf{z') is found by replacing e, 
and (7 by e^,/, rmin,wf and awf in Eq. (jB3P below. 



Appendix B: The interaction between two interfaces 

For long-ranged interparticle forces the dominant con- 
tribution to the interaction between two interfaces sepa- 
rated by distance / is from the tails of the interparticle 
potentials, which decay as inverse powers of distance. Be- 
low we calculate the interaction potential uj(l) between 
the walli-gas and gas-liquid interfaces, that is the extra 
free energy per unit area of interface for a layer of gas 
thickness /, compared to the free energy for the two sepa- 
rate interfaces (Z — >■ oo). An equivalent procedure may be 
followed to calculate the interactions between the other 
pairs of interfaces. Using the sharp-kink approximation, 
in which the fluid densities are taken to be uniform and 
equal to the bulk densities pg and respectively, on 
each side of the sharp gas-liquid interface, we flnd 
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uj{l) = {pi ~ Pg) [ pyji / v^f{z)dz ~ Pg I v{z)qIz ) (Bl) 

h+dw J I 



where pgv{z') is the potential due to a semi-infinite slab of the gas a distance z' away: 



oo poo 



Pgv{z') = pg dz 2TTr(j)att{Vr^ + z^)dr. (B2) 

Jz' Jo 

Performing the integration we obtain 

""(^^^l, fr^ z' J r-L„ , 2.1= '^^'^ r ,6 1 A , (B3) 

47re — T ^ p- + ^-5- 



4 12 6 ' 9r9 . 3r3 



Similarly, Pwi^wfiz') is the potential due to a semi-infinite slab of the wall a distance z' away where Vwf{z') is found 
by replacing e, r^in and cr by e^f, rmin,wf and cr^/ in Eq. (|B3p . Integration of (|Bip gives 



(B4) 



I 

where the temperature dependent coefficients h{T) and where bp = {pi 
c{T) are given by 

h{T) = {pgca^ - Pwif-wfcrljf) bo (B6) 

and 

c(r) = 2dwp^ie„/CT^^&o (B7) 
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